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ABSTRACT 

Based  on  arbitrarily  right-censored  observations  from  a  prob¬ 
ability  density  function  f°,  the  existence  and  uniqueness  of  the 
maximum  penalized  likelihood  estimator  (MPLE)  of  f°  is  proven. 

In  particular,  the  "first  MPLE  of  Good  and  Gaskins"  of  a  density 
defined  on  [O,00)  is  shown  to  exist  and  to  be  unique  under  arbi¬ 
trary  right-censorship.  Furthermore,  the  MPLE  is  shown  to  be  in 
the  form  of  the  solution  to  a  linear  integral  equation  whose 
forcing  function  is  an  exponential  spline  with  knots  at  the  observ¬ 
ed  censored  and  uncensored  data  points. 

1.  INTRODUCTION 

The  problem  of  nonparametric  probability  density  estimation 
has  been  studied  for  many  years.  Summaries  of  results  for  com- 


piece  (uncensored)  random  samples  have  been  listed  by  Tapia  and 
Thompson  (1978),  Wertz  and  Schneider  (1979),  and  Bean  and  Tsokos 
(1980),  for  example.  Also,  a  review  of  results  for  censored 
samples  has  been  given  by  Padgett  and  McNichols  (1984) .  In  addi¬ 
tion  to  its  importance  in  theoretical  statistics,  nonparametric 
density  estimation  has  been  used  in  hazard  analysis,  life  testing, 
and  reliability,  as  well  as  in  the  areas  of  nonparametric  discrim¬ 
ination  and  high  energy  physics  (Good  and  Gaskins,  1971). 

One  approach  to  estimating  a  density  function  nonparametri- 
cally  is  that  of  maximum  likelihood.  Nonparametric  maximum  like¬ 
lihood  estimates  of  a  probability  density  function  do  not  exist  in 
general.  That  is,  the  likelihood  function  for  a  complete  sample 
is  unbounded  over  the  class  of  all  possible  densities.  However, 
by  suitably  restricting  the  class  of  densities,  a  nonparametric 
maximum  likelihood  estimator  (MLE)  may  be  found  within  the  restric¬ 
ted  class.  For  complete  samples,  the  maximum  likelihood  estimator 
of  a  density  g  was  given  by  Barlow,  Bartholomew,  Bremner  and 
Brunk  (1972)  if  g  was  assumed  to  be  either  decreasing  (nonin¬ 
creasing)  or  unlmodal  with  known  mode.  Wegman  (1970a, b)  assumed 
unimodality  with  unknown  mode  and  found  the  MLE  of  the  density  and 
studied  its  properties  for  complete  samples.  McNichols  and  Padgett 
(1982)  studied  the  nonparametric  MLE  of  monotonic  or  unlmodal  den¬ 
sities  based  on  arbitrarily  right-censored  observations.  Even  with¬ 
in  the  class  of  decreasing  (or  unlmodal)  density  functions,  however, 
when  the  largest  observation  was  censored,  McNichols  and  Padgett 
(1982)  had  to  restrict  their  estimator  to  a  finite  interval  [0,T] 
where  T  was  an  arbitrarily  large  positive  number,  greater  than 
the  largest  observation. 

Another  approach  to  the  problem  of  nonparametric  maximum 
likelihood  estimation  of  a  density  from  complete  samples  was  pro¬ 
posed  by  Good  and  Gaskins  (1971).  This  method  allowed  any  smooth 
integrable  function  on  the  interval  of  interest  (a,b)  (which  may 
be  finite  or  infinite)  as  a  possible  estimator,  but  added  a 
"penalty  function"  to  the  likelihood.  The  penalty  function  penal- 


lzed  a  density  for  its  lack  of  smoothness,  so  that  a  very  "rough" 
density  would  have  a  smaller  likelihood  than  a  "smooth"  density, 
and  hence,  would  not  be  admissible.  De  Montrlcher,  Tapia,  and 
Thompson  (1975)  shoved  that  the  natural  mathematical  setting  for 
the  solution  of  the  maximum  penalized  likelihood  eetlmation  (MPLE) 
problem  of  Good  and  Gaskins  (1971)  was  provided  by  the  Sobolev 
subspaces  of  the  Hilbert  space  l^OO*  the  square-integrable  func¬ 
tions  on  the  real  line  R.  They  proved  existence  and  uniqueness 
results  for  the  MPLE.  Later,  Klonlas  (1982)  obtained  the  strong 
consistency  of  the  MPLE  of  the  density  function  in  appropriate 
norms.  He  also  derived  the  "first  MPLE  of  Good  and  Gaskins"  for 
the  case  that  the  density  g  has  support  only  on  the  half  line, 
essentially  by  reflecting  g  around  zero  and  using  results  for 
g  having  support  R. 

In  this  paper  we  obtain  existence  and  uniqueness  results  for 
the  nonparametrlc  MPLE  of  a  density  g  based  on  arbitrarily 
right-censored  observations  from  g.  General  results  are  first 
obtained  for  densities  with  support  ft  c  R  and  penalty  function 
$  and  then  the  problem  of  "Good  and  Gaskins'  first  MPLE"  is  con¬ 
sidered  for  arbitrarily  right-censored  data  observed  on  R.  The 
existence  and  uniqueness  results  are  then  obtained  for  densities 
g  with  only  positive  support  by  using  a  syimaetry  argument,  re¬ 
flecting  g  about  zero,  and  then  utilizing  the  general  results 
for  support  R.  It  is  also  shown  that  the  MPLE  is  in  the  form  of 
a  solution  to  a  linear  integral  equation  whose  forcing  function  is 
an  exponential  spline  with  knots  at  the  data  points. 

2.  NOTATION  AND  BASIC  DEFINITIONS 

Let  ft  c  r  be  a  finite  or  infinite  interval  and  let  f° 

denote  a  probability  density  function  with  support  in  ft.  Let 

X?,...,X°  be  n  independent  identically  distributed  random 
a  n  o  o 

variables  with  coanon  density  f  .  Later,  X^,  1*1,..., n,  will 

be  interpreted  as  the  true  survival  times  of  n  items  or  indi- 


viduals  under  observation,  where  f  will  have  support  in  [0,“). 

Suppose  that  U,  ,U_,...,U  is  a  sequence  of  constants  or  random 
12  n 

variables  which  "censor"  X^,  i*l,...,n,  on  the  right.  In  sur¬ 
vival  analysis  or  reliability  studies,  the  U^'s  represent  pos¬ 
sible  "loss"  tines  of  items  or  individuals  from  the  test. 

The  observed  data  are  denoted  by  the  pairs  (X^.A^), 
n,  where 


mlnCx' 


w-  *,  ■  j ! 


if 


x°  S  u 
i  i 


if 

,o 


>  UJ 


i  "i 
based  on  these  observa- 


It  is  desired  to  obtain  the  MPLE  of 
tlons. 

o 

In  reliability  or  survival  analysis,  where  f  has  support 
in  [0,®),  the  nature  of  the  censoring  depends  on  the  "  ' 
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CD  if  u. 


,U  are  fixed  constants,  the  observations  are  time 
’  n 


truncated.  If  all 
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are  equal  to  the  same  constant,  then  the 


case  of  Type  I  censoring  results. 

rth  order  statistic  of  X° . X° 

—  l  n 

Type  II  censoring,  (iii)  If  U 
sample  from  a  distribution  H 
ent  of  X^,...,X°,  then 


(ii)  If  all  -  X“r),  the 
then  the  situation  is  that  of 
.  .,U  constitute  a  random 


1 - n 

(usually  unknown)  and  are  lndepend 


(JC^,A^),  i»l, . . . ,n. 


is  called  a  ran¬ 
domly  censored  sample.  See  Gill  (1980,  Ch.  3  and  Ex.  4.1.1)  for 
further  discussion.  An  observed  value  of  (X^,A^)  will  be 
denoted  by  (x^d^). 

By  Lp(ft)  we  will  mean  the  space  of  functions  v  such  that 
/fl|v  (t)lpdt  <  «  with  norm  |lv||p  -  [/^I v(t)  I Pdt ]  1>^P  for  p  2  1. 
Let  H(ft)  be  a  manifold  in  L^(ft). 

Following  notation  similar  to  that  of  De  Montricher,  Tapia, 
and  Thompson  (1975),  let  ♦  denote  a  functional  $:  H(H)  •*  R. 
Given  the  arbitrarily  right-censored  sample 
the  ♦-penalised  likelihood  of  v  c  H(ft)  is  defined  by 


(xi»di) *  • • • *®» 


1-d. 


n  d  _  _ 

L(v)  -  II  [v(xJ)J  1  fl-V(x.)J  1  exp(-4»(v)). 


where  V(x.)  ■  /  v(t)dt  denotes  the  cumulative  distribution 
function  with  density  v  and  $  is  the  penalty  function.  By  the 
maximum  penalized  likelihood  estimator  (MPLE)  of  f°  correspond¬ 
ing  to  manifold  H(ft)  and  penalty  function  $,  ve  will  mean  any 
solution  to  the  problem: 

maximize  L(v)  subject  to  (2.1) 

v  e  H(n),  v(t)dt  ■  1,  and 
v(t)  2  0  for  all  t  eft. 

The  function  L(v)  is  the  censored  form  of  the  penalized  likeli¬ 
hood  of  Good  and  Gaskins  (1971) . 

When  H(ft)  is  a  Hilbert  space,  a  natural  penalty  function 
to  use  is  $(v)  *  | ! v 1 1  ,  where  ||»||  is  the  norm  on  H(f}).  If 
no  reference  is  given  to  $  when  ve  are  considering  the  MPLE 
corresponding  to  a  Hilbert  space  H(H),  it  is  assumed  that  <f>  is 
the  square  of  the  norm  on  H(f)) .  A  Hilbert  space  inner  product 
will  be  denoted  by  <*,*>  so  that  <v,v>  *  ! |v| |  .  When  H(H) 
is  a  Hilbert  space,  it  is  a  reproducing  kernel  Hilbert  space  (RKHS) 
if  point  evaluation  is  a  continuous  operation,  that  is,  yn  *  v 
in  H(fl)  Implies  that  v^(t)  -*•  v(t)  for  all  t  e  fi.  See  Goffman 
and  Pedrick  (1965)  for  further  details. 


3.  EXISTENCE  AND  UNIQUENESS  OF  AN  MPLE 

In  this  section  we  establish  the  existence  and  uniqueness  of 
a  solution  to  problem  (2.1)  when  H(ft)  is  a  RKHS.  The  inner 
product  on  H(fl)  is  defined  by  <u,v>  *  / ^  u(t)v(t)dt  for 
u,v  e  H(&). 

Theorem  3.1.  Assume  that  H(ft)  is  a  RKHS,  integration  over  ft 
is  a  continuous  functional,  and  D  is  a  closed  convex  subset  of 
{v  c  H(ft) :  v(xi)  2  0,  i»l,...,n)  with  the  property  that  D  con¬ 
tains  at  least  one  function  which  is  positive  at  the  data  points 
x1(...,x  .  Then  the  MPLE  of  f°  corresponding  to  penalty 
function  $(v)  •  |  |v|  ]  in  (2.1)  exists  in  D  and  is  unique, 
where  ||*}|  denotes  the  norm  on  H(ft). 
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The  proof  of  Theorem  3.1  is  omitted.  It  is  analogous  to  the 

proof  of  Proposition  2.1  of  De  Montricher,  Tapia,  and  Thompson 

.  -  n 

(1975),  using  the  inequality  L(v)  £  ||v||  exp(-||v|j  )(  II  K  ), 

i-1  1 

where  k  -  d^  is  the  number  of  uncensored  observations  and 

is  such  that  |v(x^)|  £  K^| |v| {  for  each  i«l,2,...,n.  Also, 
the  first  and  second  Frlchet  derivatives  of  J(v)  *  In  L(v)  are 
given  by  (Tapia,  1971) 


n  d.  r,(x.) 


(1-0  Xi 


J'(v)C-i)  -  Z  \  ±  -  I  JZT^T)  I  -.(t)dt  -  2  <  v,p.  > 

i-1  ( i  i-1  1  k{ V 

n  V*u.»  .  a'd1)I^r(t)dt|2 


J"(v)(n.n) 


n  d  n  (x  )  n 

■  Z  -K - —  -  I  - 

i-1  v  (x.)  i-1 


[l-VCx^]' 


-  2  <n,->. 


We  note  that  the  constraints  in  (2.1)  define  a  closed  convex 
subset  of  (v  e  H(H):  v(x^)  i  0,  i-l,...,n).  Also,  let  (a,b)  be 
a  finite  interval.  For  each  integer  s  i  1,  let  H*(a,b)  denote 
the  Sobolev  space  of* functions  on  Ia,b)  whose  s-1  derivatives 

are  absolutely  continuous  and  vanish  at  a  and  b  and  whose  s  th 

2  s 

derivative  is  in  L  (a, b) .  The  inner  product  on  H  (a,b)  is 

o 

defined  by 


<u,v> 


ft  u(s)(t)v(s)(t)dt. 


where  uV8‘/  denotes  the  sth  derivative.  It  is  well  known  that 

H*(a,b)  is  a  RKHS  with  the  above  inner  product  and  integration 
o 

over  (a,b)  is  a  continuous  operation  (Lemma  2.1  of  De 
Montricher,  Tapia,  and  Thompson,  1975). 

Corollary  3.1.  The  MPLE  corresponding  to  H*(a,b)  with  $(v)  - 
<v,v>  -  I [ v | | ^  exists  and  is  unique. 

As  a  special  case  of  Corollary  3.1,  we  can  consider  the  MPLE 
of  a  lifetime  density  f°  over  a  finite  interval  [0,T]  for 
very  large  T  >  0  based  on  an  arbitrarily  right-censored  sample 
from  f°.  The  MPLE  exists  and  is  unique  in  H*(0,T)  with 
penalty  function  $(v)  -  /  [vv  '(t)]  dt.  The  extension  to 


.  v  .1 


*  a  *  O  «  "  »  '  i  *  *  •  ‘  '  •  '  »  ’  •  *  •  '  *  .  * 


(0,®)  Is  considered  In  the  next  section. 

4.  THE  FIRST  ESTIMATOR  OF  GOOD  AND  GASKINS  UNDER  CENSORING 


For  complete  samples.  Good  and  Gaskins  (1971)  considered  the 
penalty  function 


♦W  ■  «  J1  ot. 


for  a  >  0,  which  is  equivalent  to 


♦  M  '  43  /l  2  <"• 


De  Montricher,  Tapia,  and  Thompson  (1975)  indicated  that  the  under 

lying  manifold  for  the  MPLE  with  this  penalty  function  should  be 
L  i  1 

v  e  H  (-“,*),  where  H  (-“,*)  is  the  Sobolev  space  of  functions 

f:  R  •*  R  such  that  the  first  derivative  f'  exists  almost  every- 
2 

where  and  f,f'  c  L  (-“,*)  with  inner  product 

<f.g>  -  f(t)g(t)dt  +  f,(t)g,(t)dt. 

% 

Letting  u  *  v  ,  we  have  the  penalty  function 

<J>(u2)  *  4a  [u'(t)]2dt,  u  e  H1(-*,«). 

This  substitution  avoids  the  nonnegativity  constraint  in  problem 

(2.1). 

For  the  data  (x^d^,  1*1,..., n,  described  in  Section  2,  we 
now  would  like  to  maximize 

L(u)  -  IT  [u2(x.)]  i[T  u2(t)dt]  1  exp(-4a|  |u’  | |2). 
i-1  1 

Since  L(u)  2  0,  maximizing  L(u)  is  equivalent  to  maximizing 
*  \ 

L(u)  *  [L(u)J  .  Thus,  we  have  the  problem: 

~  n  d„  ,  (l~d.)/2  2 

Maximize  L(u)  *  IT  [u(x  )]  1  [j  u  (t)dt]  exp(-2a] |u* | | _) 

i-1  1  xi 


subject  to  u2(t)dt  -  1. 

A 

Letting  J(u)  *  In  L(u),  problem  (4.1)  is  equivalent  to: 


(4.1) 
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Maximize  J(u)  *  I  d  In  u(x  ) 
i-1  1  1 

n  1  r*i  2  r“  2 

+  E  ±<l-d.)ln[l-/J^u  -  2^/^ [u ' (t)  j  dt 

i-1  1 

subject  to  ]_" u  (t)dt  -  1.  (4.2) 

Theorem  4,1.  Problem  (4.2)  has  a  unique  solution  in  the  set 

S  -  (u  «  H1  (-“>,«):  r  u2(t)dt  -  l). 

Proof :  The  first  part  of  the  proof  is  similar  to  arguments 
in  the  proof  of  Proposition  3.3  of  De  Montricher,  Tapia,  and 
Thompson  (1975),  but  the  details  are  somewhat  different.  The 
Frechet  derivatives  of  J(u)  are 

n  d  r,(x  )  n  (1-d.)  Xi 

j,(u)(T1)  ,  z  -L— L-  -  I  _i_  /  u(t)rt(t)dt 

i=l  uvA^'  1*1  iy  — 00 

CD 

-  4a  /  u'(t)n,(t)dt, 

-DO 

00 

2 

where  U^x^  "  /  u  (c)dt,  and 
Xi 

n  d  n2(x  )  n  (1-d  )  *i  , 

J"(u)  (n ,n)  -  -  I  -  E  - (U,(x  )  /  rf  (t)dt 

i*l  u  (x^)  i*l  (l^C^))  -00 


+  2 rf  u(t)r(t)dt)}  -4a  /  r,’(t)Zdt. 

-OO  -*X> 

Since  J"(u)  is  negative  definite,  J  is  strictly  concave,  and 
by  Theorem  2,  page  160,  of  Tapia  and  Thompson  (1978),  J(u)  has 
at  most  one  maximizer  in  the  set 

S'  -  {u  e  H1(-«,®):  u2(t)dt  s  l}. 

If  J(u)  is  continuous  on  S',  by  Theorem  4  on  page  162  of 
Tapia  and  Thompson  (1978),  J  will  have  at  least  one  maximizer 
in  S'.  To  show  this  continuity,  we  note  that  by  properties  of 
a  RKHS,  if  u  ■»  u  as  m  ®  in  (-“,»),  then 
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u  (x. )  ■*  u(x. )  for  each  i*l,  Also,  |ju  —  u  |  |  -►  0  as 

n  1  i  n  i 

m  -*■  ®  implies,  by  definition  of  the  norm  in  "  (-®,®),  that 

I  |u  —  u  |  |  _  -*■  0  and  ||u*  -u'|L-*-0  as  m  -*■  ®.  Furthermore, 

m  2  ®  -  j  ^  fa  2 

for  any  fixed  constant  c,  u^(t)dt  u  (t)dt  as  m  -*■  ®. 

Hence,  J:  S’  -*■  R  is  continuous.  Therefore,  J(u)  has  a  unique 

maximizer  u.  in  S'. 

foo  2  2 

Next,  suppose  that  u^(t)dt  <  1.  Since  u^(t)  -*■  0  as 

t  •+  ®,  then  u*(t)  and  u*(t)  both  converge  to  zero  as  t  -*■  ®. 
Thus,  there  exists  a  number  M  such  that  u^(t)  <  1  for  t  >  M. 
Consider  a  function  vA(t)  defined  so  that  (i)  v  (t)  =  u^(t) 
for  t  <  M,  (ii)  v*(t)  >  u^(t)  for  t  >  M  and  t)dt  *  1, 

and  (iii)  [v^(t)]2  <  [u^(t)]2  for  t  >  M.  Then  by  (i)  and 
(iii)  /^0tv^(t)]2dt  =  /!Llu;(t)]2dt  and 

.  Also,  by  (ii),  for  each  x^ 

1-1,..., n,  f°  v2(t)dt  >  f°  u2(t)dt.  These  results  imply  that 
xi  xi 

J(uA)  <  J(v^),  a  contradiction,  since  u^  is  the  unique  maximizer 
of  J  in  S'.  Therefore  u2(t)dt  *  1,  completing  the  proof. 

Ill 

Now,  we  assume  that  f°  is  a  lifetime  density  on  the  half¬ 
line  R+  ■  (0,®)  and  use  a  symmetry  argument  about  zero  to  obtain 
the  results  for  f°.  Thus,  assume  that  the  censored  sample 
(Xi »Af) ,  i-l,...,n,  is  such  that  >  0  with  probability  one. 
Then  the  problem  (4.1)  becomes: 

-  n  d  x  ,  %(l-d) 

Maximize  L(u)  =  II  [u(x  )]  1[l-/_(^  u  ( t ) d t ] 
i-1  1 

*  exp[-2a  /^(u'(t))2dt) ,  (4.3) 

where  xt  >  0,  1-1,..., n,  /^u2(t)dt  -  1,  and  u(t)  i  0,  t  >  0. 

Let  X_i  *  Xi  and  d_i  -  d^  1-1,..., n,  and  define 

u(x)  -  u(|x|)  for  x  e  R\{0)  and  u(0)  *  lim  u(x).  Then 

x-*0 

define  the  following  problem: 
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/.  r 


.  s. / 


Maximize  L(u)  *  II  (u(x  ))  [2 -/  *  u  (t)dt] 

|i|-l  1 

xexp[-2a  /^,(u'(t))2dt],  (A. A) 

where  u2(t)dt  -  2  and  u  e  Hg  =  {g  £  H1  (-»,»):  g(x)  *g(-x)} 

*  2 

Notice  that  L(u)  •  (L(u)3  .  Also,  H  is  equivalent  to  the 

1  b 
Sobolev  space  H  (0,®). 

Proposition  A. 2.  If  u*  solves  (A. A),  then  u*(t)  ■  u*(t), 
t  2  0,  and  u*(t)  =0,  t  <  0,  solves  (A. 3). 

~  -  u 

Proof:  Suppose  u*  solves  (A. A).  Since  L(u)  *  [L(u)K 
and  u*  is  symmetric  about  zero  implies  that  /^[u*(t)]2dt  *  1, 
u*  solves  (A. 3).  Ill 

From  Proposition  A. 2,  the  "first  MPLE  of  Good  and  Gaskins" 

2 

under  arbitrary  right-censorship  will  be  given  by  (u*)  (t).  We 
next  show  that  this  solution  exists  and  is  unique. 

Theorem  A. 3.  Problem  (A. 3)  has  a  unique  solution. 

Proof:  H„  defines  a  closed  convex  subset  of  H1^-®,®). 

Thus,  by  a  proof  similar  to  that  of  Theorem  4.1,  problem  (A. A) 

has  a  unique  solution.  By  Proposition  4.2,  u*  is  the  unique 

solution  to  problem  (A. 3).  Ill 

To  discover  the  general  form  of  the  unique  solution  u*  of 

problem  (4.4),  we  consider  the  following  problem: 

For  given  A  >  0  and  a  in  (4.1),  let 

<J>^(u).=  23jr^aj[u ' (t)  ]  2dt  +  X/^eou2(t)dt. 

n  d,  x.  Js(l-d  ) 

Maximize  L^(u)  “  II  [u(xi)3  [2 -J_^  u  (t)ot) 

li  1-1 

x  exp[-4>^(u) },  (4.5) 

subject  to  u  e  H  and  f°  u*"(t)dt  *  2. 

w 

The  inner  product  <u,v>  ■  Zaf^u'  (t)v’  (t)dt  +  u(t)v(t)dt 
defines  a  norm  ( | u | |2  ■  4>^(u)  equivalent  to  the  original  norm 
on  H*(-®,®).  Let  v,  denote  the  representer  in  the  -inner 
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product  of  the  continuous  linear  functional  given  by  point  evalu¬ 
ation  at  x^,  that  is  "  ^(x^  for  all  r>  e  h1 (-»,»). 

Let  S  ■  {v  c  Hg!  v(x^)  2  0} .  Then  S  is  closed  and  convex 
Letting  ■  In  L^,  we  have  the  first  and  second  Frechet 
derivatives, 


n  d,  r)(x. )  n  (1-dJ  *i 


E  -1—  -L-  -  Z  J  u(t)ri(t)dt  -  2<u,n>  . 

(i|-l  u(xi)  |i[-lU2(xi;  -»  X 


2 

where  U^x^  =  /  u  (t)dt,  and 


J^'(u)(n,n) 


n  d  T)  (x.)  n  U-a  ;  _ 

-  z  — —  -  z  - ^  fu  (x  >  /  n  (t)dt 

Uhi  u2(Xi)  Uhi  nyx^i  1  -> 


(1-d .) 


+  2  CJ  13(t)n(t)dt]2}-2<r!,n>^. 

Thus  - J is  uniformly  positive  definite  relative  to  S.  This 
implies  that  -J^  is  uniformly  convex  on  S.  Therefore,  if  we 
can  show  that  is  continuous  on  S,  by  Theorem  6,  page  162 

of  Tapia  and  Thompson  (1978),  will  have  a  unique  maximizer 

in  S.  By  an  argument  similar  to  that  in  the  proof  of  Theorem  4.1 
is  continuous  on  S,  and  has  a  unique  maximizer  u^  in  S. 
Now,  at  the  solution  u^,  we  must  have  the  gradient  of 
vanish.  Let  g.  be  the  element  of  H^(-®°,®)  such  that 

fx‘-  -  i,’, 

»  J  u.  (t)r}(t)dt.  Then  V  J.  (u.)  •  Z  - 

-«  i  «1  u,(xj 


n  (l-d^ 


X'"i' 
-2, 


.  ?  ~ it — F7T  “  2  u.  *  0,  where  l’.  (x  )  *  /  u.  (t)dt.  Bence 

|i|«l  2XUU  A  ZA  1  x±  X 


1 

°A  *2 


8, 


LUI-1  ux(x|)  "  ! if-!  (1'dl) 


(4.6) 


In  order  to  obtain  the  form  of  v^  in  (4.6),  from  the 


12 


definition  of  the  X-inner  product,  we  have 

2a r  v'(t)n'(t)dt  +  xf  v  (t)n(t)dt  *  n(x  ).  (4.7) 

4  —cc  i  “ 00  1  1 

Integrating  the  left-hand  aide  of  (4.7)  by  parts  (in  the  distri¬ 
bution  sense)  gives 

f ^  [XVi(t)  -  2av^(t)]n(t)dt  -  J"^  6t(t)n(t)dt  ,  (4.8) 

where  6  (t)  ■  6  (t-x  )  and  6  denotes  the  Dirac  delta  function. 

1  o  i  o 

Equation  (4.8)  is  equivalent  to  the  differential  equation 

Xvt(t)  -  2av»(t)  «  <5.(t)  (4.9) 

which,  for  i*0,  has  the  solution 


vQ(t)  -  (2aX)'%  exP[-(-^)%lt|],  t  *  0. 

Hence,  v^t)  *  v^t-x^  +  v^t+x^  solves  (4.9). 

Next,  to  determine  the  form  of  gj^  in  (4.6),  replacing  the 

x. 

right-hand  side  of  (4.7)  with  JJ^  ux(t)ri(t)dt  yields  the  non- 
homogeneous  differential  equation 

g^’(t)  -  (X/2a)gl(t)  -  -  (2a)"1  ^(t)!^  j(t),  (4.10) 

where  IA(t)  denotes  the  indicator  function  of  the  set  A.  Using 
the  theory  of  Green’s  functions,  the  solution  to  equation  (4,10) 
is 


8l(t)  -  exp  [-  tj  +  exp  [(yj)'4  t] 

-  (SaX)*’4  lC0  »lnh  [<^)!l(t-T)]l(_iiiii(!Ti);x(T)dT 


Substitution  of  vi(t)  and  g  (t)  into  equation  (4.6)  gives 
the  result  that  u^  solves  the  linear  integral  equation 


U)  .  c(ti  x.e.i)  +  (SaX)'*  /'  [  I  1 C— .^1 ( 1 T ! J] 


Ul' 

t(X/2a)%(t-T)3  Hx(T)dT, 


x  sinh 


(4.11) 


where  the  forcing  function  is  defined  by 


(  n  d  A2a\)~'*  . 

C(t;  x,a,X)  :  E  -  [exp(-(X/2aK  |t-x 

2  Mil-1  uA(xi)  1 


D 


+  exp(-(X/2a)5*  Jt+xJ)] 


-  E  .7--7-A-  [exp(-(X/2a)t)  +  exp((\/2a)\)  ] 

UI-1  U2XVV 

Note  that  over  the  constraints  of  problem  (4.4),  problems 
(4.4)  and  (4.5)  have  the  same  solution  for  any  X  >  0  since 
f*  u2(t)dt  is  constant.  By  an  argument  similar  to  that  in  Tapia 
and  Thompson  (1978,  p.  113),  the  unique  solution  u  to  (4.4) 
also  must  solve  (4.5)  with  a  positive  X,  and  thus,  has  the  same 
form  as  the  solution  to  (4<11).  The  unique  solution  to  problem 
(4.3)  is  then  u*(t)  ■  u*(t),  t  >  0,  from  (4.11),  and  hence, 
the  "first  MPLE  of  Good  and  Gaskins"  is  (u+)2. 


5.  CONCLUSION 

In  this  paper  we  have  shown  the  existence  and  uniqueness  of 
the  MPLE  of  a  density  function  in  an  appropriate  general  mathe¬ 
matical  setting,  based  on  arbitrarily  right-censored  observations 
from  that  density.  For  the  first  penalty  function  of  Good  and 
Gaskins  (1971),  the  existence  and  uniqueness  of  the  MPLE  of  the 
density  function  on  (0,»)  was  also  shown  for  this  type  of  data. 
This  "first  MPLE  of  Good  and  Gaskins"  under  arbitrary  right-cen¬ 
soring  was  shown  to  be  in  the  form  of  a  solution  to  a  linear 
integral  equation.  These  results  are  analogous  to  the  complete 
sample  case,  except  that  the  form  of  the  penalized  likelihood, 
and  therefore,  the  MPLE,  is  complicated  by  the  terms  involving 
the  survival  function. 

Statistical  properties  of  the  MPLE  under  censoring  have  not 
been  considered  here.  The  consistency,  which  seems  to  be  quite 
difficult  to  prove,  and  other  statistical  properties  need  to  be 
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investigated.  The  problem  of  computation  of  the  KPLE  from  the 
general  form  (4,11),  at  least  approximately.  Is  currently  under 
study  by  the  authors. 
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